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}_( ' Abstract 

The deformation of a fluid-fluid interface due to the thermocapillary stress induced by a continuous 
Gaussian laser wave is investigated analytically. We show that the direction of deformation of the 
liquid interface strongly depends on the viscosities and the thicknesses of the involved liquid layers. 
We flrst investigate the case of an interface separating two different liquid layers while a second part 
is dedicated to a thin film squeezed by two external layers of same liquid. These results are predictive 
f"|~| . for applications fields where localized thermocapillary stresses are used to produce flows or to deform 

O ' interfaces in presence of conflnement, such as optofluidics. 
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When a Gaussian laser beam heats the interface separating two superimposed optically absorbing liquids, 

pH I the inhomogeneous heating induces a local variation of the interfacial tension which in turn generates 

Q ■ tangential shear stresses and sets up hydrodynamic flows known as thermocapillary or Marangoni flows 

O I [Ij. These laser-induced thermocapillary flows were studied theoretically, experimentally and numerically 

in various cases, with a particular emphasis to the free-surface configuration. In the eighties, Loulergue 

et al. [2J investigated theoretically the heating of a flat fluid interface by a spatial sine modulation of an 

^ ■ infrared laser beam to produce an infrared image converter. At the same time, Viznyuk et al. [3J, inves- 

^^ I tigated the free-surface deformation of thin layers by Gaussian laser beams and applied their theoretical 

results to beam shaping using what they called optocapillarity. Beyond optical applications, Longtin et 

al. [3] investigated time dependent behaviors of free-surfaces by numerical simulation and intended a 

^^ - comparison between predictions for laser pulse heating with experiments and scaling analyses, the goal 

^— V I being to get new insights in laser melting and welding [5J. One can also cite the theoretical analysis of 

r\i ■ Rivas [1] who studied the effect of the viscous stress on the interface deformation and the numerical work 

of Marchuk [7] who investigated the effect of convection. 

We know that the temperature variation of the surface tension is usually negative for classical liquids [5]. 
Therefore, fluids are almost always pushed toward coldest regions and the previous experiments showed 
that laser heating produces beam centered dimple at the interface. This led Bezuglyi et al. [9J to consider 
^ I very thin layers and investigate the possibility to optically control fllm rupture and hole formation in or- 

" " der to understand new mechanisms for varnish dewetting, wetting [TU] and spreading [TT] of liquid drops 

and microfilms over solid substrates. However, contrary to usual expectations, Misev [121 experimentally 
observed a curvature inversion of the interface deformation, i.e. a concave to convex transition, when 
decreasing the thickness of the fluid layer. Therefore, it appears that dimple formation does not represent 
a general rule when fluid confinement starts to play a role, as already supported by classical Marangoni 
experiments in thin cells [13^. This transition from concave to convex interface deformation has been 
studied recently [14] but the generalization to two-liquid systems still deserves to be investigated. This 
is the goal of the present work. 

We theoretically investigate thermocapillary effects induced by a Gaussian laser beam that locally heats 
the interface separating two liquid layers or the interfaces of a thin film bounded by two liquid layers. Be- 
yond the general description of flow patterns and interface deformations in many different situations, our 
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investigation also provides new insights for digital optofluidics applications [15], [16], where light- actuated 
droplets are naturally confined and often close together, squeezing by the way thin liquid layers. 
Our study is structured as follows: the theoretical resolution of the flow and the interface deformation 
is detailed in section 2. Section 3 reports results and discussion on the dependence of the flow and the 
deformation to the viscosity and layer thicknesses ratios in two different configurations. We first consider 
an interface separating two different liquid layers and then generalize these results to the case of thin film 
separating two same liquid layers. 



2 Two-Fluid theoretical model 
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Figure 1: Schematic representation of the liquid-liquid system heated by a laser beam. See text for other 
notations, ct^ indicated that the interfacial tension increases with temperature. In this configuration, 
we arbitrarily consider da/dT > 0, a situation often encountered in microfluidics due to the presence of 
large amount of surfactant. 

Let us consider two liquid layers separated by an interface initially flat, horizontal and at rest, crossed 
by a continuous Gaussian laser beam propagating perpendicularly to the interface (see figure [T]). If at 
least one of the layers absorbs light, the interface will be locally heated by the laser. Physical properties 
of the liquids (denoted 1 for the bottom liquid and 2 for the top one) are their viscosities rji, 772 and 
densities pi, p2- The interfacial tension is denoted by <j. Liquids are enclosed in a cylindrical cell of 
radius R » w^ where w is the radius of the temperature distribution due to laser heating. 
Considering the axisymmetry, along the z axis of the exciting beam, of the temperature distribution, 
cylindrical coordinates (er, 69,62) with origin O located at the intersection of the beam axis with the ini- 
tial flat interface are chosen for this study. A point x is thus referenced by the space coordinates (r, 9, z). 
Two different configurations are investigated. The first one (a) considers two liquid layers separated by 
an interface, and the second one (b) generalizes the two-layer situation to a symmetric system composed 
of three liquid layers, liquid 1 being bounded on top and bottom by liquid 2 (see figure [TJ. 



2.1 Heat equations 

The temperature distribution within the fiuid layers due to the heating sources T^ obeys the heat equation 



dT 

--^ + (Ui • v)r, =. Dt, V^T,; + T, z = 1, 2 
ot 

Ui and Dxi being the velocity and the thermal diffusivity of fluid i. 

Assuming, without the loss of generality, that only fluid 1 absorbs light, T2 = and 

Ti = -—^ail{r,z) 



(1) 



(2) 



where Ai and ai are respectively the thermal conductivity and the optical absorption of fluid i at the 

used optical wavelength. 

Light intensity distribution is given by 

. ^ 2P /-2r2 \ 

I(r, z) — — TT exp — g aiz (3) 

where P is the beam power and ujq the radial extension of the laser beam, also called beam waist, cjq 
is always smaller than w due to the non locality and the strong dependence in boundary conditions of 
temperature distributions. 

The keypoint for producing thermocapillary flows is not the overheating itself but the amplitude of the 
temperature gradient. Laser heating is thus very appealing because it is extremely easy to produced weak 
amplitude overheating with large gradients. Thus we can confidently consider a weak optical absorption 
in fluid f such as aiz ^ aiHi -C 1, and neglect the axial variation of light intensity and therefore 

/W-^expf^) (4) 

The convective and unsteady terms of equation ([l} can be neglected when considering that the ther- 
mal Peclet number is small {Pct — — — ^ 1, where uq is a characteristic velocity) and a characteristic 

heating time uP' / Dti small compared to the viscous characteristic time 'uP' /vi, where Vi is the kinematic 

viscosity of fluid i. 

An example of experiment fitting these assumptions can be found in a recent publication |17) . 

Exact solutions of the heat equation using continuity of temperature and heat flux at the interface can 

be calculated P7]. However for the sake of simplicity we will consider a Gaussian distribution for Ti{r) 

which is a reasonable approximation of the exact temperature field. 



2.2 Fluid equations of motion 

Both fluids obey the mass conservation and Stokes equations, and are coupled through stress balance in 
addition to the continuity of velocity at the interface described by its height h{r) and denoted Si. As 
flows are expected to be viscous, the boundary value problem can be expressed as 

V.Ui = ; i== f,2 (5) 

= -Vp, + 77,Aui ; 1 = 1,2 (6) 

(Ti-n-T2-n) = (aK(r)-(pi-p2).g2)n+— t ; x e 5/ (7) 

ui = U2 ; X e Si (8) 

Pi is the pressure in each fluid and Ti = —pil + rii{\7ui +* Vuj) is the stress tensor corrected with a 
gravity term. The unit vector normal to the interface directed from fluid 1 to fluid 2 is denoted by n, 
s is the arc length on the interface, and t is the unit vector tangent to the interface such as (t,e6),n) is 
orthonormal. 

f d r— 
In Equation JTl), kM ~ — ; '' is the double mean curvature of the axisymmetric interface in 

da . 
cylindrical coordinates. — — is the tangential stress due to a spatial variation of the interfacial tension. 

OS 

When this variation is due to a radially inhomogeneous heating, it is usually called thermocapillary stress. 
The variation of interfacial tension can be induced by the local heating and by the migration of surface 
active surfactants molecules at the interface. Generally, when there is no surfactants, da/dT < and 
interfacial tension is smaller at the hot spot but when the interface is charged with surfactants, the effec- 
tive change of interfacial tension, involving coupling between temperature and concentration, can make 
interfacial tension larger at the hot spot. In the present investigation, we consider the effective change of 
interfacial tension da/dT due to both the heating and the solutal effect. More details on the change of 



behavior of interfacial tension with temperature, due to the presence of surfactants, can be found in the 
following references: [T7], \TE\ . 

The thermocapillary stress induces local flows on both sides of the interface which lead to its deformation. 
Considering on the one hand small deformation amplitude of height h(r) (i.e. dh/dr << 1), we can write 
the tangential stress condition as 

dui r I du2 r , 8(7 da da dT 

where T{r) is the temperature distribution at the interface due to laser heating. On the other hand, 
assuming a Gaussian axisymmetric temperature distribution such as 

T(r) = AToe^'^'/'"', (10) 

This expression of T{r) allows us to solve the hydrodynamics problem using the Fourier-Bessel transform 
defined such as 

/•OO 

h{r) ^ / h''{k)Jo{kr)kdk (11) 

where h'' is the Fourier-Bessel transform of h{r), k the reciprocal variable, and Jo{x) the 0-order Bessel 
J function. We deduce 

/•CO 

T{r) = I T^{k)Jo{kr)kdk, (12) 

2 2i.2 



with 

w 



T\k) = ATo^e-^^ (13) 

As we have — ^^-^ = —kJi{kr) where Ji{x) is the 1-order Bessel J function, we finally find 

^^ = - r T''ik)Ji{kr)k^dk (14) 

or Jo 

Therefore, we find from equation ^ that Ui^rir, z) necessarily takes the form 



Ui^rir,z)~ / Uj^ ^{z)Ji{kr)kdk (15) 

Moreover as the mass conservation condition (equation ([5])) can be written as 

i%^ + ^ = (16) 

r or oz 

it yields 

/>oo 

u,^z{r,z)= u^^^{z).h{kr)kdk (17) 



JO 

and the following relation between axial and radial velocities 

du'^ Jz) , 

-^ = -kulM (18) 

From Equations psp . (fTTl) and (fTHl) we can now calculate the velocity field due to localized laser heating. 

2.3 Velocity field 

As far as one considers small interface deformation amplitude, the velocity field can be evaluated for the 

initially flat interface h{r) = 0. 

The stokes equations (equation ([6])) can be written as follows in each liquid layer i — 1,2 

_^+^Y^ + iA_l + ^^y^^(^,2)^0 (19) 

dr \ dr^ r dr r^ dz"^ J 

- 1 -'• (1^ -;l -a »•■'<-) ^» <-) 



Remembering the following properties of the Bessel functions 

cfJi{kr) \d.h{kr) 1 



+ '-^-^Ji{kr) = -eJi{kr) (21) 



dr'^ r dr 



r 



dr^ r dr 

we rewrite Equations ([TO)l and (PUJ) as 



d^Mkr) , ^dMkr)^_^,j^^^^^ ^22) 



^+r ""'{&- ^') <rir)Mkr)kdk = 



(23) 

(24) 



Eliminating the pressure and using equation (jlSp . we finally get 



^-fc') «tW-0 (25) 

This 4*'' order linear equation of order 4 defines the eigen modes of the axial velocity. 
Its general solution is of the form 

u'^^iz) = A.e'^' + B.e-'"' + Q^e'=^ + A^e-'^^ (26) 

tii tli 

Using equation (fT8|) , the solution for the radial velocity is written 

u1,{z) = -^.e''- + B,e-'=- _ ^(fcz + l)e'=- - lTr^~kz + 1)^6"*=^ (27) 

ktii kiii 

where Ai, Bi, Ci and Di are constants to be determined from boundary conditions. We thus have eight 
unknowns to integrally solve the velocity field. The eight boundary conditions are: 
Vanishing axial velocities at the interface (z = 0) 

ulJz^0) = 4jz = 0)^0 B.C.I, 2 (28) 

Continuity of the radial velocity at the interface (z = 0) 

ul^{z ^ 0) = ui^{z = 0) B.C.?, (29) 

Vanishing velocities at the top z — H2 and bottom z = —Hi of the container (except for the symmetric 

dz 

u'l,{z = -Hi)=0 B.CA (30) 



du 

configuration (b), where it is — j-^ which is null a.t z = —Hi 



dii 

u'l^{z^-Hi)^0 B.C.ha or -^{z^-Hi) = B.C.hh (31) 

dz 

ul,{z = H2) = B.C.6 (32) 

M^,,(z = i/2) =0 B.C.7 (33) 

and finally the tangential stress jump at the interface (equation ^) for /i ~ 

dui,r du2j- da dT 

Vi^^^z=o-V2-?^\z=o^-;^jr- B.C.8 (34) 

oz oz 01 or 



B.C.I and 2 yield: 



Ai = -Bi , A2 = -B2 B.C.I, 2 (35) 



We can thus simplify notations and keep only 6 unknowns {Ai, A2, Ci, C2, Di,D2). 
B.C. 3 yield: 



-.4|-^ir--^4t4t-" 



B.C. 4 and 6 yield 



Ai{e-^"' - e''"') ~ Cie-''"' - Die''"' = B.CA (37) 

A2{e''"^ - e'''"^) + Cze*-'-^^ + Dje"''-^' = B.C.6 (38) 

B.C. 5 yield 

- Ai(e-'=^i +e'=^i) - \^{-kHi + l)e-'=^i - \^{kHi + l)e''"' = B.C.^a (39) 

k Hi k Hi 

or 

Alkie-''"' - e''"') + §-{-kHi + 2)e-''"' - ^{kHi + 2)e''"' ^ B.CM (40) 

Hi Hi 



and B.C. 7 leads to 



1<^2,,„ , .,kH2 1-02. ,„ , .^-kH2 



A^ie""' + e'""') - j^{kH2 + l)e''"- - --^{-kH^ + l)e-''"' = B.C.I (41) 

k H2 k H2 



Finally B.C. 8 leads to 



2^^{C2-D2)-2^^{Ci-Di) = -^kTHk) B.C.8 (42) 



Boundary conditions 1 — 8 entirely determine the velocity field produced by a localized laser heating and 
allow to calculate the resulting pressure field. Note that when Hi = H2, we have Ur{~z) — Ur{z) and 
Uz{^z) = —Uz{z) whatever ?7i/?72- 

2.4 Pressure field 

The pressure gradient along the z-axis (equation fiOp ) can be written as 



^ = y°° ??. (^ - kA ulMMkr)kdk 



(43) 



Integration allows to determine the pressure in fluid 1 and 2 by specifying a reference pressure PH2 at 

Z = H2 



P2{r, z) = p (pH2 + 2;^(C2e'=" + Z^ze"'^^)') Mkr)kdk 



(44) 



and 

Pi(r, z)^ j (pH2 + 2-^(Cie'=^ + Die-''')\ Mkr)kdk (45) 

which finally leads to 

P2{r,0) ~ piir,0) = J (2^^(02 +D2)-2^^{Ci+Di))Mkr)kdk (46) 

2.5 Interface deflection 

Knowing velocities and pressure fields, the shape of the interface can be deduced from the normal stress 
equation (equation ^) which can be re- written as 

'Pi{r,z) + 2r]i^^) -[-P2{r,z) + 2rn^^) ^ aArh - (pi ~ p2)gh (47) 



Considering as previously the small deformation amplitude case, the left hand side of equation (|47|) can 
be approximated by its value at z = 0. It yields 

4(7yiAi - ri2A2)kJo{kr)kdk ^ aArh - (pi - p2)gh (48) 

Using the definition of the Fourier-Bessel transform (Equation (fTT]) )and the relation A^/i = — Jq k^ h'^ (k) jQ{kr)kdk, 
where A^ denotes the radial part of the Laplacian operator, Equation (|47l) becomes 

4(7/2^2 - mAi)k = (afc2 + (p, - p2)g)h\k) (49) 

which leads to the expression of the induced deformation of the interface 

= 4 r ^-^^^^p^Mkr)k^dk (50) 

Jo Crfc + (Pl - P2)g 

Defining the gravitational Bond number as Bo = — — p^'^"^ ^ we finally find 

„, f^ A2{k)-^AAk) 



h{r)^Aw^ „,„ "\ Jo{kr)wk'dk (51) 

In the following, we use the heating length w '^ 10 — 100 fim as a reference length, uq = — i^l^jATo as 
a reference velocity and we define the thermal parameter a such as 

da ATo ,.„, 



dT a 

dT 



Assuming ??, - 1 - lOOmPa.s, a - lO^^^^-i^ ^y^ _ lOA' and |f - lO-'^Nm-^R-^, we find a - 1 
and uq ~ 5 — 5007TT.m/s. For all the calculations, we set Bo — 0.05. 
In the next section we present the velocity field and the interface deformation predicted analytically for 
two liquid layers separated by a flat interface (configuration (a)). We investigate the influence of the layer 
thickness to the heating spot ratio Hi/w, the layer ratio H1/H2 and the viscosity ratio 772/'?i on both 
thermocapillary flows and interface deformation. 

3 Results and discussion 

3.1 Double layer system 

3.1.1 Flow pattern 

An example of velocity field produced by the laser heating at the interface of a double layer configuration 
is reported in figure [5] for the case a > 1 . We observe the development of two contrarotative toroidal 
eddies induced by the tangential stress at the interface and bounded vertically by the interface and the 
walls. 

The reduced radial and axial velocity, respectively Ur/uQ and Uz/uq, are reported in figure[2]as a func- 
tion of z. We can notice that the maximum radial velocity is located at the interface {z — 0) and that the 
axial velocity is larger in the bottom layer which is thicker. Figure [3] also shows that the deformation of 
the interface is directed towards the thicker layer when fluid viscosities are the same. In the case a < 0, 
flows and interface deformation are reverted. 



3.1.2 Influence of the heating length 

Figure U shows that the maximum radial and axial velocity are increasing functions of Hi/w for 772 = ?7i. 
Note that when Hi — H2 and 771 = 772, there is no deformation of the interface because the resulting 
thermocapillary stress is null. We observe that when Hi << w, both components depend on Hi/w, as 
the radial velocity Ur scales like uqHi/w and the axial velocity u^ scales like uq{Hi/w)'^ . However, when 
Hi » w both radial and axial components saturate and scale like uq, the characteristic thermocapillary 
velocity. 
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Figure 2: Steady flow pattern in a double layer configuration for Hi = 5H2 = 5w, a = 1 and 772 — r]!- 
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Figure 3: Reduced radial velocity Ut/uq (Top left), reduced axial velocity Uz/uq (Top right) and interface 
deformation (Bottom) for Hi == 5H2 = 5w, a = 1 and 772 = Vi- 
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Figure 4: Variations of the radial (continuous) and the axial (dashed with symbols) velocities versus the 
fluid layer reduced thicknesses Hi/w {Hi ~ H2 and r — Q.lw). The right inset shows the reduced location 
zq/w of null radial velocity and maximum axial velocity (zq = z{ur = 0) = z{uzmax))- a — \ and ry2 = "Hi- 
The left inset shows the reduced deformation height h/w as a function of Hi/w for 772 — lO^yi- 



The characteristic lengths of the flow pattern are reported in the right inset of figure SI We notice that 
the largest amplitude of the axial velocity always coincides with a null radial velocity. The location of this 
characteristic point of the fiow zq scales like Hi when Hi << w as the liquid is bounded by a wall while 
it scales like w when w << Hi as the flow magnitude significantly decreases beyond the heating spot size 
w. The left inset of figure H] shows the variation of the interface deformation amplitude as a function of 
Hi/w when the liquid viscosities are contrasted (772 = lOf^i). We can first notice that when a > and 
Hi ~ H2 the deformation is always directed towards the less viscous fiuid. Moreover, the deformation 
amplitude decreases when increasing the layers thickness Hi so that Hi >> w leads to /i — ?> 0. Even 
though not strictly quantitative we can approximate this decrease so that {hl/w '^ w/Hi when Hi << w. 



3.1.3 Influence of the liquid layer thickness ratio 

Figure [5] shows the variation of the maximum axial velocity ratio Uzi/uz2 as a function of Hi/H2. We 
chose Hi = w to compare one of the heights to the length scale of the thermal forcing. We first observe 
a saturation of Uz-^/uz2 when Hi » H2 or H2 » Hi and that the axial velocity is always larger in 
the thicker layer. The inset of figure [5] shows that when 771 — r/2, the deformation is directed towards the 
thicker layer. When Hi << H2 (i.e. the heating characteristic length w = Hi is much smaller than H2); 
h only depends on w, while in the opposite case, when w » H2, the reduced deformation amplitude 
h/w strongly depends on H1/H2 = w/ H2. At the same time, the deformation amplitude h increases as 
H2 decreases as previously shown in figure |H 




Figure 5: Ratio of maximum axial velocities in each fluid Uzi/uz2 as a function of H1/H2 in a double layer 
system. Inset shows the reduced deformation height h/Hi as a function of Hi/H2- Hi = w, r = O.lw 
and r]2 = f]!- 



3.1.4 Influence of the viscosity ratio 

An interesting feature of the thermocapillary flow in our configuration is the symmetrical properties of 
the velocities when Hi = i?2- In that case, we have Ur{—z) = Ur{z) and Uz{—z) = —Uz{z) for any value 
of rj2/rji. Figure El shows that u^ax ^ uq and therefore Umax ^ I wl^2o/(2 < 77 >). This results can be 
retrieved using equation Q and scaling arguments for Hi » w. In this case —duri/dz ~ dur2/dz '^ 
Umax/w and —da/ds ^ 4^ATq/w. Similarly, using scaling arguments and equation (j47p . we find that 
Arh '^ h/w^ ^ 2(?7i — 'r]2)u„iax / (wcr) and therefore h/w ~ 0(771 — 772)/ < 7/ > as illustrated in the inset 
of figure ini 
Many features of the flows and interface deformation induced by the laser heating when da/dT > 
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Figure 6: Reduced maximum velocity Umax/uo as a function of viscosity ratio 772/771 with ryi = cste. The 
inset shows the reduced deformation height h/w as a function of 1/(1 -|- 112/ rii). Hi = i?2- 
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can be concluded from this section. First, the flow magnitude increases with the layer thickness to the 

heating length ratio Hi/w when Hi « w before reaching a saturation when Hi » w. When viscosities 

are equal in both layers, the velocity is larger in the thicker one and the interface deformation is always 

directed towards this same layer. 

Moreover, when layer thicknesses are equal. Hi = H2, the deformation is always directed towards the less 

viscous fluid. Its amplitude increases when decreasing the layers thicknesses as the vertical confinement 

increases the efficiency of the localized thermocapillary forcing. Finally, keeping Hi — H2, the dependence 

of the velocity magnitude and interface deformation can be predicted using simple scaling arguments that 

yield respectively to u^ax ^ \§^\ATq/{2 < 77 >) and h/w - |r^(?7i - m)/ < ?? >• 

In the next section, we investigate the laser heating of two fluid interface in a triple layer symmetric 

configuration where a thin film of liquid is bounded by two layers of a different liquid. 

3.2 Triple layer system 
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Figure 7: Steady flow pattern in a triple layer configuration for Hi = 0.2w, H2 = w, a — 1 and 772 — tji. 

The fiow pattern showing the generation of eddies, in a system consisting of a thin film bounded by 
two external liquid layers of same thickness, is reported in flgure [7] We can observe the production of 
two toroidal contrarotative eddies in the central film in addition to a toroidal eddy in each of the external 
liquid layers. The continuity of the radial velocity makes each pair of adjoining eddies contrarotative. 
In this section, we investigate the influence of the fllm thickness to the external layer thickness ratio 
H1/H2 on the produced thermocapillary flows and interface deformations. 

3.2.1 Influence of the film thickness 

The maximum radial velocity reduced by its value at the interface as a function of H1/H2 is reported in 
flgure [S] for equal viscosities (771 =772)- We notice that the maximal radial velocity is constant for thin 
fllms [Hi << H2) as Urmax — 0.5ur{z = 0) whilc this velocity decreases when the film thickness becomes 
comparable to the external layer thickness. 

The inset of figure [S] shows that the deformation is directed towards the external layers when Hi << H2 
and its amplitude h depends on Hi while its direction changes for H1/H2 ~ 0.7. 

It is interesting to study this system in both cases where the viscous tangential stress is directed towards 
the heating spot {da/dr < 0) or in the opposite direction {da/dr > 0). In the first case, the deformation 
of the interface is directed towards the external layers inducing a dimple (figure ^ while in the second 
case, the deformation is directed towards the thin film inducing two noses (flgure llOp . 
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Figure 8: Reduced maximum velocity Urimax / Ur{z = 0) as a function of H1/H2 in a triple layer system. 
Inset shows the reduced deformation height h/Hi as a function of Hi/H2. Hi = w and 772 = f^i- 
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Figure 9: Reduced radial velocity Ur/uo (Top left), reduced axial velocity Uz/uo (Top right) and interface 
deformation (Bottom) in a three-layer system for Hi = O.lw, H2 = 2w, a = 0.1 and 772 ~ Vi- 



3.2.2 Inducing a dimple 

When the effective interfacial tension increases with temperature, the tangential stress is directed towards 
the heating spot. Assuming liquid layers of equal viscosities 771 = r/2, the deformation of the interface 
reported in figure [5] is directed towards the external liquid layers of larger thicknesses. This prediction 
has been observed experimentally in a system similar to our configuration. In a recent investigation. 
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Dixit et al. [19] used an infrared laser beam to induce the coalescence between two SDS-coated water 
drops in Decanol. The laser beam heats the rear of one of the drops and thermocapillary stresses force 
its migration towards the other drop. The laser heating produces an accumulation of surfactant at the 
rear of the drop [2^ and therefore a deficit at the front which is in contact with the second non-heated 
drop. This deficit increases interfacial tension and therefore creates tangential stresses directed towards 
the front of the drop near the second drop. Coalescence would then occur at dimple edges location. Even 
though not strictly comparable, we believe that the dimple observed in this experiment is similar to that 
predicted by our analytical model and can be explained using the same arguments. The coalescence of 
liquid drops has also been observed experimentally in a different configuration [2T] where the heating laser 
intercepts the front and rear interfaces of water drops in contact, flowing in oil inside a microchannel, 
however the mechanism of the coalescence was not investigated in details. We believe that a dimple is 
formed in the thin film separating the drops before coalescence, similarly to the observation of Dixit et 
al. [Tg. 
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Figure 10: Reduced radial velocity UtJuq (Top left), reduced axial velocity Uz/uq (Top right) and interface 
deformation (Bottom) in a three-layer system for Hi = O.lw, H2 = 2w, a = —0.075 and 772 = '/i- 



3.2.3 Inducing noses 

In the case where the effective interfacial tension decreases with temperature, the tangential stress is 
directed from the hot region to the non-heated area of the interface. The direction of the flow is inverted 
and the deformations are directed towards the thin film inducing a nose (or hump) at each interface (see 
figure fTU|). This could also lead to drop coalescence by forming a liquid bridge similar to the observation 
made in the experimental investigation led by Bremond et al. |22) even though the mechanism at work 
here is quite different. 
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4 Conclusion 

We investigated the flows and the deformation of a liquid-liquid interface induced by the heating of a 
Gaussian laser beam in confined multi-layer systems. A first part was dedicated to an interface separating 
two liquid layers in the case where the effective interfacial tension increases with temperature. When the 
eflfective interfacial tension decreases with temperature, the direction of velocity and interface deformation 
is simply reverted. The analytical resolution showed that the flow magnitude increases when increasing 
the thickness of the liquid layers before reaching a saturation when the thickness is much larger than the 
heating lengths scale. In the case of equal layer thicknesses, the deformation is always directed towards the 
less viscous fluid and its amplitude increases when decreasing the layers thicknesses. When the viscosities 
of the layers are equal, the velocity is always larger in the thicker layer and the deformation is always 
directed towards this same layer. In the special case where both viscosities and thicknesses are equal, 
there is no deformation of the interface. A second part was dedicated to a thin film separating two layers 
of the same liquid. We showed that decreasing the ratio of film thickness to the external layer thickness 
increases the deformation which is directed towards the external layer as far as the film is much thinner 
than the external layers. Moreover, depending on the variation of interfacial tension with temperature, 
we showed that a dimple or a couple of opposite noses can be formed in the thin film, giving some insights 
for the explanation of drop coalescence observed in recent experiments. Although simple, this analytical 
resolution is predictive and therefore advances an interesting tool to explain and predict the features 
taking place in optofluidics experiments in which fluids are naturally confined. 
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